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ABSTRACT 

We present follow-up observations with the Sunyaev-Zel’dovich Array (SZA) of optically-confirmed galaxy 
clusters found in the equatorial survey region of the Atacama Cosmology Telescope (ACT): ACT-CL J0022- 
0036, ACT-CL J2051+0057, and ACT-CL J2337+0016. ACT-CL J0022-0036 is a newly-discovered, massive 
(- 10 15 M q ), high-redshift (z = 0.81) cluster revealed by ACT through the Sunyaev-Zel’dovich effect (SZE). 
Deep, targeted observations with the SZA allow us to probe a broader range of cluster spatial scales, better 
disentangle cluster decrements from radio point source emission, and derive more robust integrated SZE flux 
and mass estimates than we can with ACT data alone. For the two clusters we detect with the SZA we compute 
integrated SZE signal and derive masses from the SZA data only. ACT-CL J2337+0016, also known as Abell 
2631, has archival Chandra data that allow an additional X-ray-based mass estimate. Optical richness is also 
used to estimate cluster masses and shows good agreement with the SZE and X-ray-based estimates. Based on 
the point sources detected by the SZA in these three cluster fields and an extrapolation to ACT’s frequency, we 
estimate that point sources could be contaminating the SZE decrement at the < 20% level for some fraction of 
clusters. 

Subject headings: X-rays: galaxies: clusters: general - cosmology: observations - cosmic microwave back- 
ground 
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1. INTRODUCTION 

The Sunyaev-Zel’dovich effect (SZE) is a small (typically 
< 1 mK) distortion of the cosmic microwave background 
(CMB) spectrum caused by the inverse Compton scattering 
of CMB photons off energetic electrons of the hot intracluster 
medium of galaxy clusters (Zeldovich & Sunyaev 1969; Sun- 
yaev & Zel’dovich 1970, 1972; for reviews see, Sunyaev & 
Zeldovich 1980; Rephaeli 1995; Birkinshaw 1999; Carlstrom 
et al. 2002). The redshift independence of the SZE makes it 
a potentially powerful tool with which to search for galaxy 
clusters, especially in the distant universe. Abundances of 
clusters probe the growth of structure and have placed use- 
ful constraints on the fluctuation amplitude, erg, and the mat- 
ter density, O-m (e.g., Henry & Arnaud 1991; Viana & Liddle 
1996; Bahcall et al. 1997; Eke et al. 1998; Borgani et al. 2001; 
Reiprich & Bohringer 2002; Schuecker et al. 2003; Henry 
2004; Mantz et al. 2008; Vikhlinin et al. 2009; Mantz et al. 
2010; Rozo et al. 2010; Vanderlinde et al. 2010; Sehgal et al. 
2011). The evolution of cluster abundance with redshift is one 
of the few probes of the growth of structure and has the po- 
tential to tightly constrain cosmological parameters and pro- 
vide insight into the equation of state of the dark energy (e.g., 
Bartlett & Silk 1994; Holder et al. 2000; Haiman et al. 2001; 
Majumdar & Mohr 2004). The biggest challenge to realiz- 
ing the cosmological potential of cluster surveys is relating 
cluster mass to an observable such as integrated SZE signal. 
Well-determined masses for even a subsample of clusters can 
significantly improve parameter constraints (e.g., Majumdar 
& Mohr 2004). 

The last few years have seen significant advances in surveys 
using the SZE. The Atacama Cosmology Telescope (ACT; 
Fowler et al. 2007), the South Pole Telescope (SPT; Carl- 
strom et al. 2011), and the Planck satellite (Planck Collab- 
oration 2011a) are, for the first time, producing catalogs 
of galaxy clusters discovered through the SZE (Vanderlinde 
et al. 2010; Marriage et al. 2010; Planck Collaboration 201 lb; 
Williamson et al. 2011). The fast mapping speeds and im- 
proved sensitivities of these instruments have enabled them to 
survey large regions of the sky with sufficient depth to detect 
massive clusters, but their limited angular resolutions (> 1') 
do not allow detailed studies of cluster astrophysics with their 
data alone. 

We have performed initial follow-up observations with 
the Sunyaev-Zel’dovich Array (SZA) of three optically- 
confirmed ACT clusters in ACT’s equatorial strip. Deep, tar- 
geted SZA observations of ACT clusters provide higher sen- 
sitivity over a broader range of cluster spatial scales than that 
of ACT. In addition, the spatial filtering of the interferometer 
provides a method of cleanly disentangling radio point source 
emission from the cluster SZE signal. ACT observational 
details and target selection are presented in Section 2. SZA 
and Chandra observations and data reduction are described in 
Sections 3.1 and 3.3, respectively. The analysis method and 
results are reported in Section 4, the implications in Section 5, 
and conclusions in Section 6. 

Throughout this paper, all uncertainties are reported at 68% 
confidence and we adopt a flat, A-dominated cosmology with 
Cl M = 0.3, Da = 0.7, and Hq = 70 km s _1 Mpc _I consistent 
with recent Wilkinson Microwave Anisotropy Probe (WMAP) 
results (Komatsu et al. 201 1, 2009). 


2. ACT TARGET SELECTION 


The Atacama Cosmology Telescope (ACT) is a 6-m diame- 
ter telescope located at an elevation of 5200 m in the Atacama 
desert of Chile. Three 1024-element transition-edge-sensing 
bolometer arrays operate at 148, 218, and 277 GHz and sur- 
vey large regions of the sky mainly in two regions, a south- 
ern strip centered at declination -52.5° and in an equatorial 
strip that encompasses the Sloan Digital Sky Survey (SDSS) 
Stripe 82 (hereafter S82, Abazajian et al. 2009). For the ACT 
instrument, observation, reduction, and calibration details see 
Fowler et al. (2007), Swetz et al. (2011), Das et al. (2011), 
and Hajian et al. (2010). Initial results from the ACT south- 
ern strip include CMB power spectra (Fowler et al. 2010; Das 
et al. 201 1), compact source (Marriage et al. 201 1) and cluster 
catalogs (Marriage et al. 2010), cluster follow-up (Menanteau 
et al. 2010), and analysis of the cosmological implications of 
CMB power spectra (Dunkley et al. 2010) and cluster yields 
(Sehgal et al. 2011). 

A matched filter method (Haehnelt & Tegmark 1996; Her- 
ranz et al. 2002; Melin et al. 2006) is used for cluster detec- 
tion to produce a cluster catalog (for ACT-specific details, see 
Marriage et al. 2010). We focus on the S82 region as op- 
tical data are available from which we measure cluster red- 
shifts and estimate cluster masses via richness and luminos- 
ity. We chose three preliminary high signal-to-noise (S/N~ 5) 
cluster detections from the S82 region for follow-up as a pi- 
lot study for SZA observations; these were ACT-CL J0022- 
0036, ACT-CL J2051+0057, and ACT-CL J2337+0016. We 
will refer to these clusters as ACTJ0022, ACTJ2051, and 
A2631, where ACT-CL J2337+0016 is a previously known 
Abell cluster. The upper row of panels of Figure 1 shows the 
SZE images of these three fields expressed as Compton-y (see 
Eq. 2) made from recent ACT maps. ACTJ0022 is a new high- 
redshift cluster found by ACT through the SZE. ACTJ0022 
was spectroscopically confirmed at z — 0.81 from long slit 
observations of the brightest cluster galaxy (BCG) conducted 
at the Apache Point and Gemini South Observatories (more 
details on these observations will be described in Menanteau 
et al. 2011, in prep., where we present the X-ray and optical 
properties of the S82 ACT clusters). The S82 data provide the 
spectroscopic redshift, 0.33, for ACTJ2051 and A2631 has a 
spectroscopic redshift of 0.27 (Struble & Rood 1999). New 
maps produced with the full ACT equatorial data set reveal 
ACTJ2051 to be lower significance (S/N - 4) than inferred 
from the preliminary maps. Table 1 lists the clusters, alter- 
nate names, and redshifts. 

3. DATA 

3.1. SZA Observations 

The Sunyaev-Zel’dovich Array (SZA) is an 8-element in- 
terferometer with 3.5-m diameter dishes operating at 30 and 
90 GHz. It has an 8 GHz bandwidth, and is located at an alti- 
tude of 2200 m in the Inyo Mountains of California. Our ob- 
servations use the 30 GHz system, which has a 10.'5 FWHM 
primary beam (field of view). Typical system temperatures 
are T sys - 40-50 K, including atmospheric contributions. Six 
of the antennas are placed in a compact configuration to max- 
imize cluster sensitivity, while two are deployed at longer 
baselines to provide higher resolution data. The spatial filter- 
ing of the interferometer allows one to disentangle the small- 
scale, positive point source emission from the large-scale, 
negative SZE signal at these frequencies (for an example in 
a similar context see Reese et al. 2002). The SZA has become 
part of the Combined Array for Research in Millimeter-wave 
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TABLE 1 

SZE Image Statistics 


Cluster 

Other Name 

z 

Timing” 


Low Resolution (r MV < 2 k/l) 

High Resolution (r MV > 2 k/1) 

? int b 

(hr) 

cr 

fy/Jy beam -1 

Beam 
) ("x") 

O" CMB 

CK) 

cr 

(/iJy beam -1 ) 

Beam 

("X") 

ACT-CL J0022-0036 


0.81 

pre 

15.7 

240 

94 x 107 

38 

230 

17x20 

ACT-CL J0022-0036 



post 

21.8 

220 

91 x 101 

38 

230 

16x20 

ACT-CL J205 1+0057 


0.33 

pre 

6.7 

350 

89 x 118 

52 

300 

15x22 

ACT-CL J205 1+0057 



post 

26.9 

230 

90 x 105 

38 

230 

16x20 

ACT-CL J2337+0016 

Abell 2631 

0.27 

post 

26.6 

240 

97 x 106 

37 

210 

17 x 19 


a Pre- and post-CARMA transition of the integration of SZA with CARMA. 
b Effective on-source integration time accounting for excised data. 


Astronomy (CARMA) and our data were taken both before 
and after the CARMA transition. 

The SZA observed the three ACT clusters between March 
and May and between August and September 2010. A bright 
quasar was observed every 30 minutes, for about 2 minutes 
integration time, to monitor the system gain and phase during 
each observation. Since the fluxes of quasars can be variable 
on timescales of days or months, these sources only serve as 
secondary calibrators. We discuss absolute calibration in Sec- 
tion 3.2. 

Pre-CARMA data are reduced using a suite of MATLAB 
routines developed by the SZA collaboration (e.g., Muchovej 
et al. 2007). Post-CARMA data are reduced using Miriad 
(Sault et al. 1995). Data are excised when one telescope 
shadows another, when cluster field data are not straddled by 
two phase calibrator observations, when there are anomalous 
changes in instrumental response between calibrator obser- 
vations, when the system temperature changes dramatically 
between integrations, or when there is spurious correlation. 

The Fourier plane is also known as the u-v plane, where 
(m, v') are the Fourier conjugate variables to right ascension 
and declination, respectively. We define the radius in the u-v 
plane as r uv — A !u 2 + v 2 , which corresponds to the baseline 
length in units of the observation wavelength, A. The con- 
figurations used for these observations produce a break in u-v 
coverage at 2 k/1, with similar sensitivities for the low (r uv < 2 
kT) and high (r uv > 2 kd) resolution data (15 short and 13 
long baselines). Table 1 summarizes the on-source integration 
times, and (for the high and low-resolution data separately) 
sensitivities and the FWHMs of the synthesized beams (effec- 
tive resolution). The rms map noise, cr CMB , in CMB tempera- 
ture units for the given synthesized beams, is also given for the 
low-resolution maps. The total integration times reported cor- 
respond to the equivalent time that the full 8-element, 8 GHz 
bandwidth SZA spent on each cluster (accounting for excised 
data). 

Data visualization is done with Difmap (Shepherd 1997). 
Identification of point sources in the field is performed with 
the high-resolution ( r m , > 2 k/l) data. Point sources are mod- 
eled to remove them from the data, and the low-resolution 
data (r MV < 2 kd) are then deconvolved to produce images of 
the cluster fields with sources subtracted. We present decon- 
volved images of the point source-removed SZA observations 
as a measure of the data quality. The model fitting is per- 
formed in the Fourier plane directly. The cluster and any point 
sources in the field are modeled simultaneously using the data 
in its entirety. Modeling of the cluster and point sources is dis- 
cussed in Section 4.1 and results are discussed in Section 4.2. 

The lower panels in Figure 1 show the deconvolved SZA 
images of the three cluster fields after removal of the point 


TABLE 2 

Pre- and post-CARMA Point Source Flux Densities 


Field 

F pre 

(mJy) 

rppost 

(mJy) 

Ratio 

ACTJ0022 

0.95 ± 0.20 

0.80 + 0.14 

1.19 + 0.30 


1.07 ±0.18 

0.65 + 0.13 

1.65 + 0.36 

ACTJ2051 

8.88 + 0.33 

5.03 + 0.17 

1.77 + 0.32 


3.66 + 1.36 

2.49 + 0.61 

1.47+ 1.05 

Weighted Average 


1.51 + 0.07 


sources. Contours are multiples of twice the rms noise in 
each map, and the color scale is in units of Compton-y. 
The FWHM of the synthesized beam (effective resolution) is 
shown in the lower left of each panel. The images all cover 
the same angular scale (1 1' on a side) and have the same color 
scale to facilitate comparison. The ACTJ0022 and A2631 
fields reveal high-significance SZA detections. The point 
sources in the ACTJ2051 field are used both to strengthen our 
calibration of the SZA data and to better assess SZE decre- 
ment contamination due to such sources. 

3.2. SZA Calibration 

Pre-CARMA SZA operations routinely monitored Mars for 
absolute calibration (e.g., Muchovej et al. 2007). In the pe- 
riod right after the SZA-CARMA merger, standard calibration 
protocols had not yet been implemented. Therefore, though 
the system gains are stable, the absolute calibration is off by 
an arbitrary multiplicative factor. We calibrate the transitional 
post-CARMA data using radio sources observed both before 
and after the integration to derive an average calibration fac- 
tor that is applied to the data. This is non-ideal because the 
sources can vary; however, we note that on average we do 
not expect the source flux densities to exhibit any particular 
trend upwards or downwards over time, so the average ratio 
of several sources pre- and post-CARMA may be expected to 
follow any changes in calibration. The cluster signal will re- 
main constant, so we can validate our calibration by fitting the 
same model to the pre- and post-CARMA data and comparing 
the central Compton-y values. 

The ACTJ0022 and ACTJ2051 fields have both pre- and 
post-CARMA data. For the calibration analysis, point source 
models are fit to the data. A cluster model is simultaneously 
fit to the ACTJ0022 data. A Markov Chain Monte Carlo 
(MCMC) analysis is performed to derive flux densities and 
uncertainties (model and fitting details are described in Sec- 
tion 4.1). Ratios of pre- to post-CARMA flux densities are 
computed. The flux densities and ratios are summarized in Ta- 
ble 2. The average ratio is 1.5 + 0.1, computed with an inverse 
variance weighting. There is a wide source-to-source varia- 
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Fig. 1. — Upper. ACT Compton-y images of the three optically confirmed clusters made from recent ACT maps. We fit (excluding the cluster region in the fit) 
and remove a quadratic polynomial from the data, smooth with a Gaussian, and resample to smaller pixels. The contours are multiples of twice the rms of each 
map (black). Also shown are contours from the match filtered map (gray) as multiples of twice the rms. Lower. Deconvolved SZA images of the three cluster 
fields in units of Compton-y. The SZA images are made with short baseline data < 2 k/l after removal of point sources. Contours are multiples of twice the rms 
of each map (for details of the SZA observations presented here, see Table 1). The color scales are in units of 1(T 5 Compton-v, and solid (dashed) lines represent 
positive (negative) signal. All figures cover the same angular scale (11' on a side) and are on the same color scale to facilitate comparison. The FWHM of the 
effective resolution (synthesized beam for the SZA) is shown in the lower left comer of each panel. For visualization purposes the SZA images incorporate only 
the low-resolution SZA data and the effective beams reflect this choice. The SZA data analysis (Section 4) uses the data in its entirety. ACTJ0022 and A2631 
both yield high S/N detections. ACTJ2051 is not detected by the SZA, which is consistent with its lower SZE signal in the current ACT maps (shown here) and 
with that expected from optical data (see Section 5). In general, the qualitative agreement between the ACT and SZA data is good. 


tion in the ratios, which range between 1 .2 and 2.0, though the 
uncertainties on the ratios are relatively large and the largest 
discrepancy is < 2 cr. Our average ratio is consistent with the 
rescaling (» 1.4) found by members of the CARMA collabo- 
ration (T. Plagge, private communication). 

The cluster SZE signal will not change, but comparisons are 
complicated by its extended structure. We also fit the same f3- 
model (fixed core radius and power law index) to the pre- and 
post-CARMA data for ACTJ0022. The ratio of the central 
Compton-y parameters also yields the same factor of 1.5. That 
the cluster normalizations yield the same calibration scaling 
as the point sources lends some support to our simple ap- 
proach. As a result, the post-CARMA data are simply scaled 
by 1.5 and used for the analysis. Note that the rms values in 
Table 1 correspond to the scaled post-CARMA data. The pre- 
CARMA data absolute flux calibration is known to better than 
10% (Muchovej et al. 2007). Given the large scatter in the 
flux density ratios, we adopt a conservative 20% calibration 
uncertainty to account for potential calibration systematics. 


3.3. Chandra Observations ofA2631 

ACTJ2337 is the Abell cluster A2631, which has publicly 
available Chandra X-ray observations. The data consist of 
two ACIS-I observations, obsIDs 3248 and 11728, of dura- 
tions 9 and 17 ks respectively. The data are reduced start- 
ing with the level 1 events file using CIAO 4.3 and calibra- 
tion database version 4.4.1. Standard corrections are applied, 
along with light curve filtering and other standard processing 
(for reduction details see Reese et al. 2010). The Chandra 
images of A2631 are made by binning the 0. 7-7.0 keV data 
to 17968 pixels and exposure maps are computed at 1 keV. 
The images and exposure maps from the two observations are 
combined, and a wavelet based source detector is used on the 
combined image and exposure map to identify and generate a 
list of potential point sources. The list is used as the basis of 
our point source mask. 

Figure 2 shows the smoothed Chandra data (color) along 
with the deconvolved SZA data (contours). The color scale 
shows Chandra counts and the contours are multiples of twice 
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Fig. 2. — Deconvolved SZA image of A2631 (contours) on the smoothed 
Chandra counts image (color). Contours are multiples of twice the rms of 
the map (see Table 1 for details) and the color scale is in units of Chandra 
counts. The FWHM of the synthesized beam (effective resolution) is shown 
in the lower left corner. The E-W elongation is seen in both the SZA and 
Chandra data, and the peaks in the SZE and X-ray data align well. 


the rms of the SZA map. The FWHM of the synthesized beam 
(effective resolution) is shown in the lower left corner. The 
same E-W elongation seen in the SZA data is seen in the 
Chandra data as well, and the peaks in the SZE and X-ray 
images align. 


4. ANALYSIS 

4.1. Method 

We fit cluster models to the data and derive parameters from 
the best-fit models. ACTJ0022 has no X-ray data and we fit 
models to the SZE data only. A2631 has both SZA and Chan- 
dra data and we perform three complementary analyses using 
1) SZA data only; 2) Chandra data only; and 3) both the SZA 
and Chandra data jointly. Derived quantities such as mass and 
integrated SZE signal are computed within Rgoo, the radius 
within which the mean interior density is 500 times the critical 
density at the redshift of the cluster, p c (z) — 3H 2 (z)/(8nG), 
where H(z) is the Hubble parameter, and G is Newton’s grav- 
itational constant. 

4.1.1. Chandra Spectroscopy for A2631 

Both a single temperature and a temperature profile are 
measured for A2631. We extract spectra and response files 
for both data sets separately. Single temperature spectra are 
extracted within (0. 15— 1 ) 7 ? 5 oo for the isothermal temperature 
analyses. The region used for spectroscopic extraction is 
found recursively, by picking a trial radius, extracting a spec- 
trum to determine the electron temperature, T e , within that 
radius, and computing M 500 = M(fGoo) from hydrostatic equi- 
librium (Eq. 3) and the resulting Rgoo- This process is re- 
peated until the values of the input and output R 500 agree to 
1%. Spectra are also extracted from within the full R 500 deter- 
mined as outlined above for comparison. 

For the temperature profile, we extract spectra in radial an- 
nuli with r out /r ln set to a constant ratio. This follows Vikhlinin 


TABLE 3 

A263 1 Chandra Spectral 
Properties 


R in 

n 

Rout 

(") 

(&) 

z 

(Zo) 

0 

309 

7 q+0.5 
' ‘ -0.4 

O 

U) 

00 

1 + 
© P 
£ £ 

46 

309 

7 q+0.6 

' -0.5 

0 42+ 013 
-0.12 

0 

20 

B.3! 2 ; 2 

1 17+0.54 

1-1 -0.55 

20 

39 

8.H!; 2 

°- 38 t26 

39 

79 

9.9!J; 2 

0.49-: 22 

79 

157 

7 0+O.8 
* ‘ -0.5 

0.48!° : 22 

157 

315 

6 7 +0 - 6 
u - 1.2 

1 77+O.36 
1 ‘ -0.36 

315 

630 

60-1:1 

0-88 


Note. — Fits performed with 
redshift z = 0.27 (Struble & 
Rood 1999) and Galactic H I col- 
umn density N H = 3.55 x 10 20 
cur 2 (Kalberla et al. 2005). 


et al. (2006, hereafter V06) who found this choice produces 
roughly equal counts per annulus for cluster observations. 
However, we choose r out /r; n = 2 (instead of 1.5 used by V06) 
in order to increase the counts in each radial bin, since our 
data have many fewer counts than the data considered in V06. 

We construct quiescent ACIS background datasets by re- 
projecting the ACIS blank-field observations to match each 
dataset. To account for variations in the particle background, 
the blank-field observations are scaled by the ratio of fluxes 
in the 9.5-12 keV band, where the Chandra effective area is 
essentially zero and the flux is due entirely to the background 
(e.g., Vikhlinin et al. 2005). Background spectra are extracted 
from these quiescent background data sets using the same re- 
gions as that for the cluster data and used when fitting the 
cluster spectral properties. 

Spectra for both Chandra observations using their respec- 
tive responses are fit simultaneously to the same plasma 
model with the normalizations allowed to vary independently 
between data sets. Xspec (Arnaud 1996; Dorman & Arnaud 
2001) is used to model the intracluster medium with an APEC 
spectrum (Smith et al. 2001) that includes Bremsstrahlung 
and line emission components. We adopt the Galactic col- 
umn density of N„ = 3.55 x 10 20 cm ~ 2 (Kalberla et al. 2005), 
solar abundances of Asplund et al. (2009), and cross sections 
of Balucinska-Church & McCammon (1992) with an updated 
He cross section (Yan et al. 1998). The analysis uses data in 
the 0.7-7. 0 keV energy range. The “cstat” statistic, which is 
similar to the Cash (1979) statistic, is used when modeling the 
spectral data to properly account for low counts. 

Table 3 summarizes the Chandra spectral results, listing the 
inner and outer extraction radii, derived temperatures, T e , and 
metallicity relative to solar, Z. In all cases, the model provides 
a good fit to the data, without any obvious structure or pattern 
in the residuals. As a check, the T e profile spectra are fit fixing 
the metallicity to the global value (Z = 0.4). The resulting 
T e ’s agree with those in Table 3 to within l<x for all regions. 

4.1.2. Cluster Models 

For the SZA data, the clusters are fit with both the tradi- 
tional /7-model (Cavaliere & Fusco-Femiano 1976, 1978) and 
the universal pressure profile of Arnaud et al. (2010, hereafter 
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A10). For the SZE-only /3-model analysis, we fix /3 = 0.86 
from a best fit to an average SZE profile (Plagge et al. 2010). 

For the Chandra data, the cluster is fit with the traditional 
isothermal /3-model and the modified /3-model (e.g., V06) for 
density along with the corresponding V06 T e (r) profile, here- 
after referred to as the V06 model. Many of the V06 model 
parameters must be fixed, given the large number of parame- 
ters compared with the quality of the X-ray data. A263 1 does 
not show signs of a cool core so we excise those parts of the 
V06 model. 

Joint fits to the SZA and Chandra data are performed by 
assuming the A10 pressure profile and a simplified, core-cut 
form of the V06 density profile (hereafter sV06) for the cluster 
models. These fits are performed both including and exclud- 
ing the X-ray spectroscopic data, allowing the SZE and X-ray 
surface brightness data to place constraints on the cluster tem- 
perature (for details see Mroczkowski et al. 2009). The tem- 
perature is derived assuming the ideal gas law ( P e = n e k B T e ). 
The X-ray surface brightness data determine the sV06 density 
fit, and the SZE surface brightness data drive the A10 pressure 
profile fit. A single temperature is derived from the inferred 
T e profile using the prescription of Mazzotta et al. (2004) over 
the (0.15-l)Rgoo region. For the A10+sV06 fits that include 
spectroscopic information, the likelihood from the inferred T e 
described above is included in the MCMC using the output 
from the Chandra single-temperature spectroscopy. A sep- 
arate MCMC run is done without including the T e spectral 
component in order to assess how X-ray surface brightness 
and SZE data can be used to find M(r) without relying on 
X-ray spectroscopy. 

4.1.3. Radio Point Source Model 

The radio point source model accounts for the primary 
beam attenuation and includes a spectral index that models 
the frequency dependence. The spectral dependence of the 
point source model is given by 



where F v is the intrinsic point source flux density at frequency 
v, F Vo is the intrinsic flux density at fiducial frequency vo, and 
a is the spectral index. We adopt the average frequency for the 
fiducial frequency, vo = 30.938 GHz, and report flux densities 
at this frequency. 

Point sources are first modeled with both F Vn and oc allowed 
to vary. The 8 GHz bandwidth of the SZA potentially pro- 
vides leverage on both the spectral index and flux density of 
each point source. The flux density is, in general, well con- 
strained. However, only weak constraints are obtained on a 
for all but the brightest sources. We first run fits including a in 
the Markov chain and then determine a using our 30 GHz flux 
density combined with 1 .4 GHz flux densities from the NRAO 
VLA Sky Survey (NVSS; Condon et al. 1998). As a check, 
we compare the flux densities obtained from the NVSS -based 
spectral indices to those obtained just from the SZA data. The 
flux densities are consistent within the 68 % confidence re- 
gions. 

4.1.4. Markov Chain Monte Carlo Analysis 

We perform a Markov chain Monte Carlo (MCMC) anal- 
ysis of the SZE and X-ray data (for details see Reese et al. 
2000; Bonamente et al. 2004). The philosophy behind the 
analysis is to keep the data in a reduced but native state, and 


to run the models through the observing strategy to compare 
directly to the data. Interferometers measure the Fourier trans- 
form of the sky brightness modulated by the primary beam. 
Therefore the SZA data most naturally provide constraints in 
the u-v plane, where the noise properties of the data and the 
spatial filtering of the interferometer are well understood. We 
perform model fits directly in the u-v plane, computing the 
Gaussian SZE likelihood, _£sze- For X-ray data, the Poisson 
likelihood, Jlx, is computed for each pixel, ignoring those that 
fall within the point source mask. The SZE and X-ray data are 
independent, and their likelihoods can simply be multiplied 
to obtain the combined likelihood. Best-fit parameters and 
confidence intervals are determined from the 50%, 16%, and 
84% levels of the cumulative distribution function, which de- 
fine the the median and 68 % confidence region. The resultant 
probability distributions are visually inspected, and conver- 
gence and mixing of the Markov chains are checked with the 
Geweke Z-statistic (Geweke 1992). 

The cluster and any detected point sources in the field are 
modeled for the SZA data. The SZE signal varies as 

AT sze = f(v)T CMB f dC-^—^n e cr T = f(v)T CMB y, (2) 
J m e c z 

where f(v) is the frequency dependence of the SZE at fre- 
quency v, n e and T e are the electron number density and tem- 
perature, k B is the Boltzmann constant, m e is the mass of the 
electron, <x T is the Thomson cross section, integration is along 
the line of sight (, and v is the Compton-y parameter. SZA 
data have sixteen 500 MHz bands covering 8 GHz of band- 
width. The frequency dependence of the SZE is taken into 
account when modeling. Relativistic corrections to /(v) (e.g., 
Itoh et al. 1998; Challinor & Lasenby 1998) depend on T e and 
are not included in this analysis for consistency, since only 
one cluster has a measured T e . The effects of this small (< 3% 
at 30 GHz) correction are discussed in Section 5. 

For A2631, Chandra data are modeled with a cluster model 
and a constant X-ray background. Regions containing point 
sources are excluded from the fit. 

4.1.5. Derived Cluster Properties 

Derived quantities such as the integrated Compton-y, T 500 , 
and mass, M500, are computed for each step in the Markov 
chain for each of the above types of fits we perform. We 
use these output chains to determine the median values and 
68 % confidence regions for each parameter of interest. This 
method cleanly propagates uncertainties from the parameters 
included in the chain. For example, uncertainties from mod- 
eling detected radio sources are propagated into the integrated 
Compton-y results. 

We compute the integrated Compton-y, K ph (r) = 
cr T /(m e c 2 ) fP e dV, within a spherical volume enclosed by ra- 
dius r, where dV - 4nr 2 dr and P e = n e k B T e is the elec- 
tron pressure. Since this latter quantity tracks thermal en- 
ergy ( E = 3/2 fP,„,dV), and thermal pressure is the domi- 
nant source of support against gravitational collapse (see, e.g. 
Nagai et al. 2007; Lau et al. 2009, who report that only 10- 
20 % of the pressure comes from non-thermal support), we 
can expect T sph (r) to track gravitational energy within radius r. 
Assuming the virial relation, a constant gas fraction of 0.13, 
an average metallicity, and a total mass profile that can be de- 
scribed by a Navaro, Frenk, and White (NFW) halo model 
(Navarro et al. 1996, 1997), we estimate mass from fits to 
K sph (r) derived from our A10 and /3-model fits to the SZE 
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TABLE 4 

Radio Point Sources 


Field 

R.A. (J2000) 
(h m s) 

Decl. (J2000) 
(d m s) 

R a 

(') 

G 1 
(mJy) 

G.4 

(mJy) 

a 

pest 
1 148 

(mJy) 

ACTJ0022 

00 22 13.006 

-00 36 33.35 

0.3 

1 oo+0.15 
^ -0.15 

18.8!^ 

-0.9 ±0.1 

0.3 


00 22 05.987 

-00 35 54.95 

1.6 

0 98 +013 
-0.13 

33.2+“ 

-1.1 ±0.1 

0.2 

ACTJ205 1 

20 51 17.758 

+00 53 20.16 

2.3 

7 36 +0 - 14 
' -0.14 

43-8Gj;t 

-0.6 ±0.1 

3.0 


20 51 40.720 

+00 52 02.74 

7.3 

4 ni + 0 - 44 
^• U1 _ 0.44 

16 9 +0J 
AO- -0.7 

-0.5 ±0.1 

1.9 

A2631 

23 37 40.209 

+00 16 42.18 

0.7 

0 qq+ 0.16 

- 0.17 

107.0+3* 

-1.1 ±0.1 

0.8 


23 37 38.106 

+00 10 01.93 

6.2 

3 20 +0 - 34 
- 0.33 

1 9 4 + 0.6 
AZ< - 0.6 

-0.4 ±0.1 

1.6 


Projected distance from the cluster center. 


TABLE 5 

SZE and X-ray Derived Properties 


Field 

/Goo 

4/500 

Goo 

(model) 

(Mpc) 

(10 14 M g ) 

(10 4 Mpc 2 ) 

actjoo22 

/?- SZE 

1 07 +004 
a,u/ -0.04 

8 -^!:o 

1.04j°; 2 ° 

A10-SZE 

A2631 

1 oo +005 

A - -0.05 

7 o+l.O 
' -1.0 

0-8411S 

13- SZE 

1 42+0.15 
1 -0.13 

1L5^;1 

1 oc+0.85 
-0.49 

A10-SZE 

1 ^o+0.36 
1 -0.25 

16.0!> 3 5 5 

1 -K-oil 

(3-X-my 

1 9 1 +0.03 
A - -0.03 

8.5j°;* 


V06-X-ray 

1 54+O.IO 
A - -0.11 

13-8+15 


A10+sV06 a 

1 90+0.20 
A,J -0.12 

Q 4+4.8 
* -2.4 

1 03 +0 ’ 31 
A - -0.17 

A10+sV06 b 

1 31+0.35 
a,ja -0.38 

O 1+9-4 
^-5.9 

1 oi +0 - 65 

1 ,u -0.52 


a Includes X-ray spectroscopic information. 
b Does not include X-ray spectroscopic informa- 
tion. 


data (as was done in Mroczkowski 2011). These SZE-only 
mass estimates can be performed regardless of chosen model 
fit, and, like estimates from the X-ray assuming hydrostatic 
equilibrium, rely on spherical symmetry and do not take into 
account sources of non-thermal pressure support. We define 
y 50 o = T sph (R5oo) as the integrated Compton-y within R$ oo- 
Mass estimates from the Chandra data are based on hydro- 
static equilibrium (e.g., Sarazin 1988), 

M( r ) - k s T e( r ) r l d Mn e ) dln(T e ) \ _ r 2 dP^(r) 
Gpm p \ d ln(r) + d ln(r) / Gp gas (r) dr 

(3) 

where p is the mean molecular weight, m p is the mass of the 
proton, and p ga , and P gas are the total gas density and pres- 
sure, respectively. The mass as a function of radius is used to 
compute R 500 using M 500 = M(R 500 ) = (47r/3)/?| 00 500p c (z), 
where p c (z) is the critical density at redshift z. Our analysis 
is similar to that of Vikhlinin et al. (2006). 

4.2. Results 

The MCMC results for all detected point sources are pre- 
sented in Table 4. In addition to the SZA 3 1 GHz flux densi- 
ties, the corresponding NVSS (Condon et al. 1998) 1.4 GHz 
flux densities are listed. These two flux densities are used 
to compute the spectral index ot, which is then used to esti- 
mate the point source flux density in ACT’s 148 GHz band, 
Fftg (both also listed in Table 4). The projected distance from 
the cluster, R, is also listed and gives an idea of the potential 
impact of each source on cluster detection and potential con- 


tamination of the SZE flux. Table 5 summarizes the cluster 
modeling results for the SZA, Chandra, and joint analyses. 
We report /Goo. M500, and T500 from our MCMC analysis. 

Figure 3 shows the SZE u-v radial profiles of the SZA 
data (points) with best- fit /? (blue) and A10 (red) models 
for ACTJ0022 (left) and A2631 (right). The visibilities ( u - 
v plane data) are converted to a frequency independent u-v 
plane Compton-y, Y uv , and scaled by the angular diameter 
distance squared, Z) 2 , creating a SZE luminosity-like quan- 
tity (for details see, e.g., Mroczkowski et al. 2009). The real 
parts of Y uv are plotted. The imaginary components are con- 
sistent with zero. Residuals are shown in the lower sections 
of both panels. The ACTJ0022 radial profile has x 2 = 35.1 
and 35.8 with 28 degrees of freedom (dof) for the / 3 and A10 
models, respectively. The corresponding probabilities of ob- 
taining these x 2, s or larger by chance given the degrees of 
freedom are, p(> x 2 1 dof) = 0.17 and 0.15. The A2631 radial 
profiles have^ 2 = 38.1 and 36.1 with 28 dof for the /? and A10 
models, respectively, with corresponding p(> x 2 \ dof) = 0.10 
and 0.14. 

Figure 4 shows the Chandra surface brightness profile (up- 
per) and temperature profile (lower) along with the best fit 
models. The data (points) and best fit / 3 (blue) and V06 (red) 
models are shown. The best-fit background level is shown by 
the dotted line. The lower portion of each panel shows the 
residuals in units of the standard deviation. Despite this clus- 
ter’s E-W elongation (Fig. 2), a spherical model, on average, 
provides a good fit with x 2 = 190.3 and 155.9 for the // and 
V06 radial profiles that have 194 and 191 degrees of freedom 
(dof), giving p(> x 2 \ dof) = 0.56 and 0.97, respectively. The 
V06 temperature profile results in \ 2 — 2.6 for 2 degrees of 
freedom with p(> x 2 \ dof) = 0.27. 

4.3. Optically Informed Cluster Properties 

SDSS S82 data provide redshift information and enable 
mass estimates based on relations between cluster mass and 
optical properties such as richness and luminosity (e.g., John- 
ston et al. 2007; Reyes et al. 2008; Rozo et al. 2009). We have 
computed cluster masses from the AGoop-AGoo scaling relation 
of Rozo et al. (2009), which is based on weak-lensing and X- 
ray measurements of clusters from the SDSS maxBCG cata- 
log (Koester et al. 2007). Masses computed from the AGoop- 
M500 scaling relation of Reyes et al. (2008) yield similar re- 
sults after correcting for the different radii used between the 
studies. M500 is the mass within /Goo. as defined earlier, the 
radius within which the mean interior density is 500 times that 
of the critical density at the redshift of the cluster. AGoop is the 
number of red sequence galaxies in a cluster measured within 
Rioop, the radius within which the mean interior density is 200 
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Fig. 3. — SZE profiles for ACTJ0022 (left) and A2631 (right) as a function of u-v radius, r uv = Vw 2 + v 2 . Best-fit A10 (red) and /5 (blue) model fits are also 
shown. Residuals are shown in the lower sections in units of the standard deviation. Point sources have been subtracted from the visibilities and the phase centers 
shifted to the cluster centers before computing the profiles. The u-v plane Compton-y, Y uv , is scaled by Z) 2 , creating a SZE luminosity-like quantity (for details 
see, e.g., Mroczkowski et al. 2009). The real parts of Y uv are shown. The imaginary parts are consistent with zero. The dotted line shows the shadowing limit of 
the 3.5 m dishes. 


Radius (kpc) 

10 100 1000 



Fig. 4. — Upper: A263 1 Chandra surface brightness profile (points) with 
best fit /? (blue) and V06 (red) models. The lower section shows residuals in 
units of the standard deviation. The models are similar and the simple spheri- 
cal /5-model provides a good fit to the Chandra data. Lower: A263 1 Chandra 
temperature profile (points) with best fit V06 temperature (red) model. The 
lower section shows residuals in units of the standard deviation. 


times the mean matter density at the redshift of the cluster, 
p(z ) = O m(z) p c (z), and is denoted by the subscript p. /Loop 
at z = 0 corresponds to /Go (60 with respect to critical), and 
is substantially larger than /Goo and /Goo- Cluster A^oop’s are 
computed for the S82 data using the maxBCG prescription as 
implemented by Menanteau et al. (2010b, see Section 2.2 for 
details). Table 6 summarizes the measured Aboop and inferred 
values of M 500 and /Goo for each cluster field. 

A recent Planck study explores the SZE-optical scaling 


TABLE 6 

Optically Informed Properties 


Field 

N 200 p 

V 500 

(x 10 14 M o ) 

^500 

(Mpc) 

I 500 

(xlO -4 Mpc 2 ) 

ACTJ0022 
ACTJ205 1 
A2631 

113 ± 11 
25 ±5 
77 ±9 

7.77 ± 1.12 
1.54 ±0.33 
5.19 ±0.75 

1.03 ±0.05 
0.73 ± 0.05 
1.11 ±0.05 

°- 84 -a!l 

Note. — 

loop’s are 

measured from SDSS S82 data, M 500 


is derived from N 2 oop using the scaling relation of Rozo et al. 
(2009), and /Goo is inferred from the mass. This /Goo is then 
used to compute Fsoo from the SZA data for comparison to the 
Planck results (Planck Collaboration 2011c). 

relations by employing a multi-frequency matched filter on 
Planck maps at the positions of the SDSS maxBCG clusters 
(Planck Collaboration 2011c). This work finds an offset be- 
tween the measured integrated C om pt on-y-op deal richness 
relation compared to model predictions for the full maxBCG 
sample. However, when using the X-ray subsample, the au- 
thors find good agreement between the prediction and the 
model. Following Planck Collaboration (2011c), we deter- 
mine .ft 500 from the optical properties via the A^oop-Afsoo rela- 
tion (Table 6 ) and compute T 500 within that radius, using the 
fits of the A10 profile to the SZA data. T 500 is rescaled to 
redshift z — 0 using the evolution of the Hubble parameter 
for a flat universe, E{z) = \JHm( \ + z) 3 + Q A , and a fiducial 
distance D A = 500 Mpc, as 

Z 500 = Y 50Q E(zy 2/3 (D A (z)/500 Mpc) 2 [arcmin 2 ]. (4) 

Figure 5 compares our two ACT clusters (red points) to those 
of Planck (black points) for the Y 500 -N 200 P relation. The best 
fit power law for the Planck data over the full maxBCG cat- 
alog is also shown (line). The known X-ray cluster, A2631, 
lies above the Planck relation, though with large uncertainty. 

5. DISCUSSION 

The density of compact radio sources is higher in cluster re- 
gions compared to the field (Cooray et al. 1998; Massardi & 
De Zotti 2004; Lin & Mohr 2007; Coble et al. 2007). Contam- 
ination of the SZE signal from radio sources could potentially 
bias flux and mass estimates from the SZE. The radio source 
fluxes extrapolated to 148 GHz (Table 4) provide information 


High-Resolution SZA Observations of ACT Clusters 


9 


on potential contamination of the SZE decrement signal. The 
two sources with projected radius from the cluster R > 6' will 
have little impact on cluster detection and flux. ACT 148 GHz 
equatorial maps have a noise level of around 2 mJy beam -1 so 
that most sources are expected to fall well below the ACT 
map noise level, and the brightest source is expected to be 1 .5 
times the noise level. Very Large Array (VLA) observations 
of galaxies in nearby clusters between 5 and 40 GHz find that 
about 60% of the point sources show a flattening of the spec- 
tral shape above 8 GHz (Lin et al. 2009). This implies that 
extrapolating from low frequency to high frequency yields a 
lower limit on the contaminating flux. A conservative upper 
limit on the radio source flux density at 148 GHz is obtained 
by increasing the extrapolated estimates by a factor of two. 
However, recent simulations of the microwave sky suggest 
that only 3% of clusters have their 148 GHz SZE decrements 
contaminated at the > 20% level (Sehgal et al. 2010). There 
is no indication in the ACT data of contamination by sources 
in these three cluster fields. The estimated 148 GHz flux den- 
sities suggest that radio sources do not significantly impact 
cluster detection in surveys for the brightest clusters, although 
they could potentially bias flux measurements in some frac- 
tion of the clusters at a level of 3-6 mJy (< 20% for typical 
-30 mJy integrated SZE fluxes). We reiterate that extrapo- 
lation to higher frequencies is uncertain and note that this is 
based on sources in only 3 cluster fields. More precise flux 
density estimates at ACT frequencies will be obtained through 
observing a larger number of ACT cluster fields over a range 
of frequencies. 

The VLA FIRST survey (Becker et al. 1995) covers the 
ACTJ0022 and A2631 fields. Three out of the 4 radio 
point sources have NVSS and FIRST flux densities that agree 
within 68% confidence. The 107 mJy NVSS source in 
A2631’s field has a 91 mJy flux density in the FIRST cat- 
alog, resulting in a predicted 148 GHz flux density that is 
9% higher than that estimated from the NVSS flux density. 
These surveys are in different configurations and observe the 
same fields at different times, thus providing a rough handle 
on potential time-variability of these sources, which could im- 
pact the contamination of the SZE signal and our use of point 
sources to calibrate our SZA data (Section 3.2). A 10% flux 
density variability is encompassed by our choice of a 20% 
calibration uncertainty for the SZA data. 

The S82 data for ACTJ2051 exhibit a red sequence typi- 
cal of galaxy clusters and provide a spectroscopic redshift. 
The optical mass estimates from the S82 data suggest that 
ACTJ2051 is less than a third of the mass of the other two 
clusters (Table 6). Using the Planck ?5oo--^200p relation, we 
estimate T500 = 2.7 + 1.2 x 1 0 6 Mpc 2 , which corresponds to 
415 ± 190 fi]y of integrated SZE flux within /Goo at 31 GHz. 
This corresponds to an SZA signal smaller than 2cr for the ob- 
servations considered here, below the SZA detection thresh- 
old. Furthermore, the above significance calculation assumes 
that the entire SZE flux within /Goo is actually contained 
within the synthesized beam (effective resolution) of the SZA. 
However, the SZA synthesized beam is T.7, smaller than /Goo 
(2'.6) for this cluster, and the cluster SZE signal will be diluted. 

Our mass estimates for ACTJ0022 and A2631 confirm that 
they are massive systems (M500 — 10 15 M Q ). This is con- 
sistent with optical and X-ray follow-up of the initial cluster 
results from ACT’s southern region that suggests the high sig- 
nificance detections will be > 8 X 10 14 M Q (Menanteau et al. 
2010). The p and A10 model fits to the SZA data yield con- 



Fig. 5. — Integrated Compton-y-optical richness relation, Y50Q-N200P, for 
Planck (black) and this work (red). Also shown is the best fit relation to the 
Planck data (solid line). Integrated Compton-}/ values are spherical values 
within R 500 scaled to redshift zero and a distance D A = 500 Mpc. Richness, 
N200p, is the number of red sequence galaxies within /?200p- 

sistent masses for both clusters. A comparison of the SZE, 
X-ray, and joint mass estimates of A263 1 shows more scat- 
ter. This scatter might be indicative of A263 1 being an unre- 
laxed cluster, especially in light of the asymmetric structure 
seen in both the SZA and Chandra observations (see Fig. 2). 
However, given the uncertainties in the masses, the different 
methodologies yield consistent results for A2631 as well. 

ACTJ0022 was discovered by ACT and we present ini- 
tial mass estimates but there are previous X-ray analyses of 
A2631. Temperatures are computed from regions that dif- 
fer between the different works. Using the earlier of the 
two Chandra observations, Maughan et al. (2008) finds T e = 
6.5 ± 0.6 keV within (0. 1 5-1 )/Goo, where /Goo = 1.2 Mpc. 
Analyses of XMM-Newton data using a number of methods 
for computing /Goo yields T e — 7.5 + 0.5 keV (within (0.2- 
0.5)/Goo) and an average M500 = (8 + 2) x 10 14 M Q within 
/Goo = 1.2 Mpc (Zhang et al. 2008). Another analysis of 
XMM-Newton data finds T e = 7.7+0.6 keV within an aperture 
containing 90% of the background-subtracted surface inten- 
sity (Andersson et al. 2009). Our derived T e = 7.3 + 0.6 keV 
agrees to within 1 .3cr with other Chandra and XMM-Newton 
measurements, despite the differences in radii. The Zhang 
et al. (2008) Mgoo is consistent with our derived SZE and X- 
ray masses to within 2cr in the most discrepant case. Despite 
data from various observatories and different methodologies, 
the temperatures and masses of A263 1 are consistent. 

A number of systematics may affect these mass measure- 
ments, summarized in Table 7 as percentages for a single 
cluster. Both simulation and analysis of X-ray data suggest 
that asphericity typically affects mass estimates at the 5-10% 
level when measured at large radii such as /Goo (e.g., Math- 
iesen et al. 1999; Piffaretti et al. 2003). We conservatively 
adopt 10% for the effects of asphericity. Simulations sug- 
gest that small-scale fluctuations, often called clumping, will 
cause a - 10% overestimate of mass (Mathiesen et al. 1999). 
Non-thermal pressure support may affect both SZE and X- 
ray mass estimates, causing mass to be underestimated at the 
10-20% level (e.g., Nagai et al. 2007; Lau et al. 2009). For 
the SZE-only analysis, changes in the assumed gas fraction, 
on average, change the mass by 12%, with lower gas fraction 
leading to higher mass; the fixed metallicity assumption is a 
< 1% effect (Mroczkowski 2011), which we neglect. Rela- 
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TABLE 7 

Potential Systematics on Mass 


Source 

SZE 

(%) 

X-ray 

(%) 

Calibration 

±10 

±10 

Asphericity 

±10 

±10 

Clumping 

-10 

-10 

Non-thermal P e 

+ 15 

+ 15 

Assumed f g 

±12 


Relativisic corrections 

+3 


Total 3 

+24 

-21 

+21 

-17 


Added in quadrature. 


tivistic corrections to the SZE are a 3% correction at 30 GHz 
for an 8 keV cluster like A2631 (Itoh et al. 1998). 

The SZA and Chandra calibrations both affect the mass es- 
timates. The SZE mass determinations depend on the SZE 
signal as M oc (A T) 1 ^ 2 so that the conservative 20% calibra- 
tion results in a 10% change in mass. Recent changes in the 
Chandra calibration can change the cluster temperatures in- 
ferred from spectroscopic fits by - 10% (Reese et al. 2010; 
Nevalainen et al. 2010). This will directly impact our X-ray 
mass estimates since M oc T e . 

Potential systematics are summarized in Table 7, with totals 
added in quadrature. Both SZE and X-ray estimates have po- 
tential systematics at the - 20% level, roughly the same order 
as the statistical uncertainties (Table 5). 

6. CONCLUSION 

We obtained SZA follow-up observations of three optically- 
confirmed clusters from preliminary ACT maps along the 
celestial equator. ACT-CL J0022-0036 is a massive, high- 
redshift cluster newly discovered by ACT. The SZA detects 
two of the three clusters at high significance. ACT and SZA 
data show good qualitative agreement (Figure 1). The clus- 
ter A2631 shows good agreement between SZE and X-ray 
data (Figure 2) with both peaks aligning and similar E-W 
elongation seen in both wavebands. Initial mass estimates 
confirm that ACTJ0022 and A263 1 are massive clusters with 
M 500 - 10 15 M q . 

Two compact radio sources are detected by the SZA in each 
cluster field (Table 4). Using NVSS 1.4 GHz flux densities, 
we compute spectral indices and predict the flux densities in 
ACT’s 148 GHz band. The radio sources are expected to be 
< 6 mJy in the ACT 148 GHz maps, suggesting that radio 
sources are not a significant contaminant for detection of high 
mass clusters. However, they can still impact the measured 
SZE signal. As an example, the brighter source in A2631 
could be filling in the SZE signal at the < 20% level, assuming 
the extrapolation of the sources’ lower-frequency flux densi- 
ties holds. A more precise determination of potential contam- 
ination of the SZE signal from compact sources at ACT fre- 
quencies will be obtained through observing a larger number 
of ACT cluster fields over a range of frequencies. 

Optimized for different purposes, ACT and SZA provide 
complementary SZE data on galaxy clusters. With the ability 
to quickly image large regions of the sky, ACT is well-suited 
to finding clusters. Targeted observations with the SZA al- 
low deep integrations for detailed cluster studies. The SZA 
images (Fig. 1) show peak S/N of 12 and 10 for ACTJ0022 
and A2631, respectively, compared to S/N - 6 for the match- 
filtered ACT images. The SZA provides higher resolution 
than ACT, measuring angular scales 15"-5', well-matched to 


clusters. In addition, the spatial filtering of the interferome- 
ter provides a clean method of disentangling cluster emission 
from radio point source emission. 

This pilot study of SZA follow-up observations of ACT- 
detected clusters shows that the detected clusters are massive 
systems. Cluster abundances are exponentially sensitive to 
mass (e.g.. Press & Schechter 1974), with the most massive 
clusters providing the most leverage on cosmological studies. 
Finding extremely massive systems at high redshift can poten- 
tially rule out the current A-CDM paradigm (e.g., Mortonson 
et al. 2011). In addition, cosmological determinations utiliz- 
ing massive systems are expected to be less susceptible to the 
effects of non-gravitational astrophysics. The highest mass 
systems from ACT’s two survey regions will comprise our 
core sample for derivation of cosmological parameters from 
clusters. It is crucial to measure well the integrated SZE sig- 
nal and the masses of these clusters for proper cosmological 
interpretation of cluster yields from current and future SZE 
cluster surveys. Deep, targeted SZA observations provide a 
method of disentangling the point source and cluster emis- 
sions and enable robust estimates of the integrated SZE signal, 
7500 , and initial mass estimates. When combined with mass 
measures from X-ray and weak lensing, this will provide a 
robust measure of the mass-SZE flux scaling relation. 


We are grateful to lohn Carpenter for guidance on CARMA 
observing, the queue, and a crash-course in Miriad. We thank 
the CARMA collaboration, especially Tom Plagge, for many 
discussions of the nitty-gritty details of the data. Support for 
TM was provided by NASA through the Einstein Fellowship 
Program, grant PF0-1 10077. The CARMA 3.5-m observa- 
tions presented here were awarded in proposals c0563 and 
c0619. 

This work was supported by the U.S. National Sci- 
ence Foundation through awards AST-0408698 for the 
ACT project, and PHY-0355328, AST-0707731 and PIRE- 
0507768. Funding was also provided by Princeton Univer- 
sity and the University of Pennsylvania. The PIRE program 
made possible exchanges between Chile, South Africa, Spain 
and the US that enabled this research program. Computations 
were performed on the GPC supercomputer at the SciNet 
HPC Consortium. SciNet is funded by: the Canada Founda- 
tion for Innovation under the auspices of Compute Canada; 
the Government of Ontario; Ontario Research Fund - Re- 
search Excellence; and the University of Toronto. 

Support for CARMA construction was derived from the 
Gordon and Betty Moore Foundation, the Kenneth T. and 
Eileen L. Norris Foundation, the James S. McDonnell Foun- 
dation, the Associates of the California Institute of Tech- 
nology, the University of Chicago, the states of California, 
Illinois, and Maryland, and the National Science Founda- 
tion. Ongoing CARMA development and operations are sup- 
ported by the National Science Foundation under a coopera- 
tive agreement, and by the CARMA partner universities. 

Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating Institu- 
tions, the National Science Foundation, the U.S. Department 
of Energy, the National Aeronautics and Space Administra- 
tion, the Japanese Monbukagakusho, the Max Planck Soci- 
ety, and the Higher Education Funding Council for England. 
The SDSS is managed by the Astrophysical Research Con- 
sortium for the Participating Institutions. The Participating 
Institutions are the American Museum of Natural History, As- 


High-Resolution SZA Observations of ACT Clusters 


11 


trophysical Institute Potsdam, University of Basel, Univer- 
sity of Cambridge, Case Western Reserve University, Uni- 
versity of Chicago, Drexel University, Fermilab, the Institute 
for Advanced Study, the Japan Participation Group, Johns 
Hopkins University, the Joint Institute for Nuclear Astro- 
physics, the Kavli Institute for Particle Astrophysics and Cos- 
mology, the Korean Scientist Group, the Chinese Academy 
of Sciences (LAMOST), Los Alamos National Laboratory, 
the Max-Planck-Institute for Astronomy (MPIA), the Max- 
Planck-Institute for Astrophysics (MPA), New Mexico State 
University, Ohio State University, University of Pittsburgh, 
University of Portsmouth, Princeton University, the United 
States Naval Observatory, and the University of Washington. 


This work made use of observations obtained with the 
Apache Point Observatory 3.5-meter telescope, which is 
owned and operated by the Astrophysical Research Consor- 
tium and observations obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under a cooperative agreement 
with the NSF on behalf of the Gemini partnership. 

This research has made use of the NASA/IPAC Extragalac- 
tic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under con- 
tract with the National Aeronautics and Space Administration. 
This research has made use of data obtained from the Chandra 
Data Archive and the CIAO software provided by the Chandra 
X-ray Center (CXC). 


REFERENCES 


Abazajian, K. N„ et al. 2009, ApJS. 182, 543 

Andersson, K., Peterson, J. R., Madeiski, G., & Goobar, A. 2009, ApJ, 696, 
1029 

Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference 
Series, Vol. 101, Astronomical Data Analysis Software and Systems V, 
ed. G. H. Jacoby & J. Barnes, 17 

Arnaud, M., Pratt, G. W., Piffaretti, R., Bohringer, H., Croston, J. H., & 
Pointecouteau, E. 2010, A&A, 517, A92 
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 
Bahcall, N. A., Fan, X., & Cen, R. 1997, ApJ, 485, L53 
Balucinska-Church, M. & McCammon, D. 1992, ApJ, 400, 699 
Bartlett, J. G. & Silk, J. 1994, ApJ, 423, 12 
Becker, R. H„ White. R. L„ & Helfand, D. J. 1995, ApJ, 450, 559 
Birkinshaw, M. 1999, Phys. Rep., 310, 97 

Bonamente, M., Joy, M. K., Carlstrom. J. E., Reese, E. D., & LaRoque, S. J. 

2004, ApJ, 614, 56 
Borgani, S., et al. 2001, ApJ, 561, 13 
Carlstrom, J. E„ et al. 2011, PASP, 123, 903 

Carlstrom, J. E„ Holder, G. P, & Reese, E. D. 2002, ARA&A, 40, 643 
Cash, W. 1979, ApJ, 228, 939 

Cavaliere, A. & Fusco-Femiano, R. 1976, A&A, 49, 137 
— . 1978, A&A, 70, 677 
Challinor, A. & Lasenby, A. 1998, ApJ, 499, 1 
Coble, K., Bonamente, M., Carlstrom, J. E., Dawson, K., Hasler, N., 
Holzapfel, W., Joy, M., La Roque, S., Marrone, D. P„ & Reese, E. D. 
2007, AJ, 134, 897 

Condon, J. J.. Cotton, W. D., Greisen, E. W., Yin, Q. F„ Perley, R. A., 
Taylor, G. B„ & Broderick, J. J. 1998, AJ, 1 15, 1693 
Cooray, A. R., Grego, L., Holzapfel, W. L., Joy, M., & Carlstrom, J. E. 

1998, AJ, 115, 1388 
Das, S., et al. 2011, ApJ, 729, 62 

Dorman, B. & Arnaud, K. A. 2001, in Astronomical Society of the Pacific 
Conference Series, Vol. 238, Astronomical Data Analysis Software and 
Systems X, ed. F. R. Harnden Jr., F. A. Primini, & H. E. Payne, 415 
Dunkley, J., et al. 2010, arXiv: 1009.0866 

Eke, V. R„ Cole, S„ Frenk, C. S„ & Patrick Henry, J. 1998, MNRAS, 298, 
1145 

Fowler, J. W„ et al. 2010, ApJ, 722, 1 148 

Fowler, J. W., et al. 2007, Appl. Opt., 46, 3444 

Geweke, J. 1992, in Bayesian Statitstics IV, ed. J. M. B. et al. 

(Oxford: Clarendon), 169 

Haehnelt, M. G. & Tegmark, M. 1996, MNRAS, 279, 545 

Haiman, Z„ Mohr, J. J., & Holder, G. P. 2001, ApJ, 553, 545 

Hajian. A., et al. 2010, arXiv: 1009.0777 

Henry, J. P. 2004, ApJ. 609, 603 

Henry, J. P. & Arnaud, K. A. 1991, ApJ. 372, 410 

Herranz, D., Sanz, J. L., Barreiro, R. B.. Hobson, M.. Martinez-Gonzalez, 

E. , & Diego, J. M. 2002, in Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, Vol. 4847, Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, ed. J.-L. Starck & 

F. D. Murtagh, 50 

Holder, G. R, Mohr, J. J., Carlstrom, J. E., Evrard, A. E., & Leitch, E. M. 
2000, ApJ, 544, 629 

Itoh, N„ Kohyama, Y„ & Nozawa, S. 1998, ApJ, 502, 7 

Johnston, D. E., Sheldon, E. S., Tasitsiomi, A., Frieman, J. A., Wechsler, 

R. H„ & McKay, T. A. 2007, ApJ, 656, 27 
Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M„ Bajaja, E., 
Morras, R.. & Poppel, W. G. L. 2005, A&A, 440, 775 
Koester, B. P, et al. 2007, ApJ, 660, 239 
Komatsu, E., et al. 2009, ApJS, 180, 330 
Komatsu, E., et al. 2011, ApJS, 192, 18 
Lau, E. T„ Kravtsov, A. V., & Nagai, D. 2009, ApJ, 705, 1 129 
Lin, Y. & Mohr. J. J. 2007, ApJS, 170, 71 

Lin, Y., Partridge. B., Pober, J. C., Bouchefry, K. E., Burke, S., Klein, J. N.. 
Coish, J. W„ & Huffenberger, K. M. 2009, ApJ, 694, 992 


Majumdar, S. & Mohr, J. J. 2004, ApJ, 613, 41 

Mantz, A., Allen, S. W„ Ebeling, H„ & Rapetti, D. 2008, MNRAS, 387, 
1179 

Mantz, A., Allen, S. W„ Rapetti, D„ & Ebeling, H. 2010, MNRAS, 406, 
1759 

Marriage, T. A., et al. 2010, arXiv: 1010. 1065 
Marriage, T. A., et al. 2011, ApJ, 731, 100 
Massardi, M. & De Zotti, G. 2004, A&A, 424, 409 
Mathiesen, B., Evrard, A. E., & Mohr. J. J. 1999, ApJ, 520, L21 
Maughan, B. J., Jones, C., Forman, W., & Van Speybroeck, L. 2008, ApJS, 
174, 117 

Mazzotta, P„ Rasia, E., Moscardini, L., & Tormen, G. 2004, MNRAS, 354, 
10 

Melin, J., Bartlett, J. G„ & Delabrouille, J. 2006, A&A, 459, 341 
Menanteau, F., et al. 2010, ApJ, 723, 1523 
Menanteau, F., et al. 2010b, ApJS, 191, 340 

Mortonson, M. J., Hu, W., & Huterer, D. 2011, Phys. Rev. D, 83, 023015 

Mroczkowski, T. 201 1, ApJ, 728, L35 

Mroczkowski, T., et al. 2009, ApJ, 694, 1034 

Muchovej, S., et al. 2007, ApJ, 663, 708 

Nagai, D., Kravtsov, A. V., & Vikhlinin, A. 2007, ApJ, 668, 1 

Navarro, J. F„ Frenk, C. S„ & White, S. D. M. 1996, ApJ, 462, 563 

— . 1997, ApJ, 490, 493 

Nevalainen, J., David, L., & Guainazzi, M. 2010, A&A, 523, A22 
Piffaretti, R., Jetzer, P., & Schindler, S. 2003, A&A, 398, 41 
Plagge, T., et al. 2010, ApJ, 716, 1 1 18 

Planck Collaboration, Ade. P A. R., et al. 2011a, arXiv: 1 101.2022 
Planck Collaboration, Ade. P A. R., et al. 2011b, arXiv: 1 101.2024 
Planck Collaboration, Aghanim, N., et al. 201 lc, arXiv: 1 101 .2027 
Press, W. H. & Schechter, P. 1974, ApJ, 187, 425 

Reese, E. D., Carlstrom, J. E., Joy, M., Mohr, J. J., Grego, L., & Holzapfel. 
W.L. 2002, ApJ, 581,53 

Reese, E. D., Kawahara, H„ Kitayama, T., Ota, N., Sasaki, S., & Suto, Y. 
2010, ApJ, 721,653 

Reese, E. D., Mohr, J. J., Carlstrom, J. E., Joy, M., Grego, L., Holder, G. P., 
Holzapfel, W. L., Hughes, J. P., Patel, S. K., & Donahue, M. 2000, ApJ, 
533,38 

Reiprich, T. H. & Bohringer, H. 2002, ApJ, 567, 716 
Rephaeli, Y. 1995, ApJ, 445, 33 

Reyes, R., Mandelbaum, R., Hirata, C., Bahcall, N., & Seljak, U. 2008, 
MNRAS, 390, 1157 
Rozo, E„ et al. 2009, ApJ, 699, 768 
Rozo, E„ et al. 2010, ApJ, 708, 645 

Sarazin, C. L. 1988, X-ray emission from clusters of galaxies (Cambridge 
Astrophysics Series, Cambridge: Cambridge University Press, 1988) 
Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, ASP Conf. Ser. 77: 
Astronomical Data Analysis Software and Systems IV, 4, 433 
Schuecker, P, Bohringer, H., Collins, C. A., & Guzzo, L. 2003, A&A, 398, 
867 

Sehgal, N., Bode, P, Das, S., Hernandez-Monteagudo, C., Huffenberger, K., 
Lin, Y.-T., Ostriker, J. R, & Trac, H. 2010, ApJ, 709, 920 
Sehgal, N., et al. 2011, ApJ. 732, 44 

Shepherd, M. C. 1997, in Astronomical Society of the Pacific Conference 
Series, Vol. 125, Astronomical Data Analysis Software and Systems VI, 
ed. G. Hunt & H. Payne, 77 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, J. C. 2001, 
ApJ, 556, L91 

Struble. M. F. & Rood, H. J. 1999, ApJS, 125, 35 
Sunyaev, R. A. & Zeldovich, I. B. 1980, ARA&A, 18, 537 
Sunyaev, R. A. & Zel’dovich, Y. B. 1970, Comments Astrophys. Space 
Phys., 2, 66 

— . 1972, Comments Astrophys. Space Phys., 4, 173 
Swetz, D. S., et al. 2011, ApJS, 194, 41 
Vanderlinde, K.. et al. 2010, ApJ, 722, 1180 
Viana, P. T. P. & Liddle, A. R. 1996, MNRAS, 281, 323 


12 


Reese et al. 


Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C., Markevitch, M., Murray, 
S. S., & Van Speybroeck, L. 2006, ApJ, 640, 691 
Vikhlinin, A., et al. 2009, ApJ, 692, 1060 

Vikhlinin, A., Markevitch, M„ Murray, S. S., Jones, C., Forman, W., 
Williamson, R„ et al. 2011, arXiv:1101.1290 


Yan, M., Sadeghpour, H. R., & Dalgarno, A. 1998, ApJ, 496, 1044 
Zeldovich, Y. B. & Sunyaev, R. A. 1969, Ap&SS, 4, 301 
Zhang, Y.-Y., Finoguenov, A., Bohringer, FI., Kneib, J.-R, Smith, G. P., 
Kneissl, R„ Okabe, N„ & Dahle, H. 2008, A&A, 482, 451 


